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The theory of cylindrical magnetic domains provides conditions governing 
the size and stability of circular cylindrical magnetic domains in plates of 
uniaxial magnetic materials together with an estimate of the range of 
applicability of these conditions. The results of the theory are directly 
applicable to the design of cylindrical domain devices. Computation to first 
and second order of the energy variation resulting from general small 
deviation in the domain shape from an initially circular shape yields the 
conditions governing domain size and stability. The physical origin of 
the various terms in the energy expansion is examined in detail. A graph 
from which many domain size and stability properties may be obtained 
summarizes the results of the energy variation calculation. The minimum 
theoretically attainable domain diameter is approximately ov/ttM 2 , , where 
<r w is the wall energy density and M, is the saturation magnetization. 
For domains to exist, the effective anisotropy field must be greater than 
4wM B . 

I. INTRODUCTION 

The recent development of a technique for the propagation of isolated 
magnetic domains in an arbitrary direction in anisotropic ferromagnetic 
thin films by P. C. Michaelis created a renewed interest in the use of 
domain propagation for device purposes. 1 The technique used by 
Michaelis for propagating domains along the easy axis is quite different 
from that used for propagation along the hard axis. During discussions 
on the possible application of these techniques, A. H. Bobeck, U. F. 
Gianola, R. C. Sherwood, and W. Shockley suggested that for general 
symmetrical domain propagation the direction of magnetization must 
lie normal to the plane of the film 2 . The recognition that rare earth 
orthoferrites have the required properties came in response to this 
suggestion. 3 Experimental work on the application of this type of 
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domain motion device was then begun. Although at the present time 
this work has been largely concentrated on the orthoferrites, there exist 
other materials, such as the hexagonal ferrites and manganese bismuth, 
having the required properties. 

The present work directs attention to structures in which the prop- 
erties of the material used require the magnetization to he normal to 
the surface of the plate. The modes of operation of devices constructed 
from such structures are classified according to the effect of wall motion 
coercivity. In the case of very high wall motion coercivity, the applica- 
tion of shaped applied fields determines the initial domain configuration 
which is then maintained by coercivity. For very low wall coercivity, 
on the other hand, the saturation magnetization, wall energy, plate 
thickness and bias field determine the domain size and shape. Between 
these two extremes, there is a continuum of intermediate modes. In 
either extremal mode, a complete set of operations (logic, memory, and 
transmission) may be performed. 4 The present work concerns only the 
low coercivity mode and specifically, right circular cylindrical domains 
in plates of uniform thickness and small variations therefrom. When 
observed by means of the Faraday effect, cylindrical domains have the 
appearance (particularly when in motion) of bubbles and therefore are 
colloquially referred to as "bubbles". 

The present work largely treats the theory of cylindrical domains with 
experiments and applications being considered only briefly. Section II 
presents the domain model and mode of description. Section III contains 
the calculation of the energy derivatives used in the investigation of 
domain size and stability. Section IV contains an interpertation of the 
energy derivatives in terms of fields and potentials. Section V discusses 
the solution of the domain size and stability equations. Section VI 
discusses the range of validity of the domain model used in the previous 
sections. It is found that several assumptions implicit in the model are 
related, and a requirement on materials suitable for the production of 
circular domains is obtained. Appendix A contains a derivation the 
properties of certain elliptic integrals appearing in the theory of circular 
domains, Appendix B is a listing of the standard forms and series 
expansions of the magnetostatic force and stability functions, and Ap- 
pendix C is a list of mathematical symbols. 

II. THE DOMAIN MODEL AND MODE OF DESCRIPTION 

Figure 1 shows the magnetic domain structure to be considered 
here. 5 The isolated magnetic domain is magnetized downward while 
the remainder of the plate is saturated upward. The domain will be 
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Fig. 1 — Magnetic domain configuration. 

considered to be near circular. Examination of the variation of domain 
energy under a variation of domain shape from the assumed unper- 
turbed shape yields domain stability. Once created, a cylindrical domain 
continues to exist if the magnetic configuration meets the conditions 
for stable equilibrium. The stability of a given configuration, however, 
does not guarantee that it can be produced. The generation of cylin- 
drical domains is a separate problem which is not treated here. 

2.1 Description oj the Domain 

A cylindrical (r, 8, z) coordinate system is placed at the center of the 
domain with its z-axis perpendicular to the plane of the plate. The 
plate is taken to have planar surfaces and a uniform thickness h. Only 
the case of a plate of infinite extent, r f = «> , is considered here. It is 
assumed that the material constraints allow the magnetization to He 
only along the z-axis and the magnitude of the magnetization is inde- 
pendent of the local magnetic field. The boundary between the two 
regions of magnetization, the domain wall, is assumed to be independ- 
ent of z (no wall bulging) and to have a width which is negligible in 
comparison to the domain radius. It is assumed that a wall energy 
density per unit area <j w may be assigned independently of either the 
orientation or curvature of the wall. The assumptions about the de- 
tailed magnetic configuration (the magnetization magnitude and orien- 
tation and the wall energy and shape) are coupled by the material 
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properties. Section VI contains a detailed discussion of the validity of 
these assumptions and the cylindrical wall assumption. Even though 
the foregoing assumptions appear quite drastic and restrictive, experi- 
mentally there does exist a region in which the results obtained under 
these assumptions are both accurate and useful. 
The expansion 

n(6) = f; r n cos [n(d - d n )] (1) 

n = 

of r b (0) in terms of the Fourier coefficients, r„ and 0„ , describes the 
domain shape in the plane. The n value is called the "rotational per- 
iodicity." The condition 



r 



» £ n | r„ | (2) 



assures that the domain is near circular and that the function r b (d) 
is single valued and smooth. 

It is convenient to introduce the finite variations of the r„ and n , 
Ar n and A0„ , respectively, in order to describe small variations in 
domain size and shape from the strictly circular domain of radius 
r [n(0) = r ]. In terms of these variations, a small variation of the wall 
shape from r b (6) = r may be written as 

00 

r b (6) = r + Ar + X) Ar n cos [n{6 - d n - A0„)] (3a) 

n=l 

where, by assumption, 

\r |» | Ar | + f>| Ar„|. (3b) 

n-l 

Subject to the restrictions stated, equation (3) describes an arbitrary 
variation because of the completeness of the Fourier expansion. 

The externally applied magnetic field, H, is taken to be spatially 
uniform and to lie in the positive z direction. (The presence of a compo- 
nent of the applied field in the plane of the plate has no effect to the 
approximation that the magnetization lies only along the z-axis.) 

The assumed simple forms of the applied field and magnetic con- 
figurations permit the use of simple formal expressions for these quan- 
tities. The expression for the externally applied field is 

h = m t (4) 

where # is a constant and i z is the unit vector in the 2-direction. The 
magnetization may be written in terms of the unit step function, 
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as 



M = Mf,{l - 2u[r b (6) - r]\u(z + \h)u(-z + \h). 
2.2 The Energy Variation 



(5) 



(6) 



The investigation of domain size and stability proceeds by computing 
the first and second variations of the total system energj' with respect 
to the r n and 0„ . The total energy of the domain is 



E t — Ejp -\- E H -f- E 



(7) 



where E w is tlie total wall energy, E„ is the interaction energy with the 
externally applied field, and E M is the internal magnetostatic energy. 
The total wall energy, under the previously stated assumptions, is the 
product of the wall energy density a w and the wall area a: 



E, 



-L 



a w da = ha w 



rl(6) 



dn(e) 
dd 



(IB. 



(8) 



The interaction energy of the magnetization with the externally applied 
field is 



E„ = - / MUdV = - 



V J -x ^0 

and the internal magnetostatic energy is 

E »-l/ F £. v ' ,MV,; "' dy ' rfV 



I\LHr dr dd dz, 



(9) 



r — r 



dz dz 7 ' 



/•;•' dr dd dz dr' dd' dr' 



where 



s a = r + r" - 2rr' cos (6 - 9') + (z - z')\ 



(10a) 



(10b) 



In expressions (9) and (10), V indicates volume and primes indicate 
quantities in the second coordinate system used in describing the in- 
ternal magnetostatic interaction. 
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The variation in the total energy when the r n and 0„ are varied is 



AE T = Z 



dE T \ . . (dEr 



?H 



I CO CO 

+ EE 

•" n = m = 



d 2 E, 



Ar„Ar, 



+ 



d 2 E, 



dd„ dd„ 



_\dr n dr„ 
A0 n AO m \ + 3 



3 2 E T 

dr n ae* 



Ar„A0„ 



(11) 



where the subscript refers to evaluation of the partial derivatives at 
the circular domain state, r 6 (0) = r , and 3 refers to terms of order 
three and higher in the combination of Ar n and A0„ . The first partial 
derivatives of the energy, (dE T /dr n ) and (dE T /dd n ) , are the generalized 
forces of the system, while the second derivatives of the total energy 
form the elements of the stiffness matrix. 

Knowledge of the generalized forces and the stiffness matrix com- 
pletely characterizes domain size and stability. It is shown in Section 
III that only the energy derivatives, (dE T /dr )o and the (d 2 E T /drl) , 
are non-zero when r 6 (0) = r . The equation obtained by setting the 
only nonzero generalized force equal to zero is called the "force; equa- 
tion." The expansion (1) is a quasi-normal mode expansion since cir- 
cular domains are completely metastable with respect to the 0„ and the 
stiffness matrix is diagonal with respect to the r n . 



III. CALCULATION OF THE ENERGY DERIVATIVES 

3.1 Derivatives of the Wall Energy 

The derivatives of the total wall energy are computed by substituting 
the wall shape expression (1) into the wall energy expression (8), 
noting that 



3d 



= — 2 nr n sin [w(0 - 0„)] 



(12) 



and differentiating under the integral sign. There results 



dE 

dr 



?-*•-/."{ 



cos [n(6 - 6„)] -j£n sin [n(0 - 0,,)] 



' [ fe 



dd 



(13a) 



and 
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d Ew 
dr„ dr„ 



= ha w 



cos [n(d - 9 n )] cos [m(6 - 9 m )] 



dn 
dd 



+ nm sin [n(0 - 0„)] sin [m(0 - 0,,.)] 
- L, cos [n{6 - 0„)] - f^nsm [n(0 - 0„)]| 
■lr h cos [m(d - 6 m )\ - ~ wi sin [m(0 - 0„)]} 



1 f*/ 



-i 



de 



(13b) 



with analogous expressions for 

dE w /dO n , d 2 E w /d6 n dd m , and d 2 E w /dr n dO, 
Evaluating equations (13) for a circular domain, 

r h (fi) = r and [dr b (0)/dd] = 0, 
the circular domain derivatives are 
f 8Bi 



dr ) Q 

d 2 E w 



= 2irh<r w 



= — h<r w n , 



n > 1 



(14a) 



(14b) 



and all of the first and second derivatives of the total wall energy not 
explicitly stated are zero. 

3.2 Derivatives of the Applied Field Interaction Energy 

The applied field interaction energy is evaluated by substituting the 
formal expressions for the applied field (4) and the magnetic configura- 
tion (6) into the applied field interaction expression (9), changing the 
order of integration, and integrating. 

/<;„ = -MM (' T [ (1 - 2u[r b (0) - r]} 

X u(z + lh)u(-z + hh)r dz dr dO (15a) 

= hM.H I rl(6)(io\ - constant. (15b) 

_* J 

The infinite constant is independent of the r„ and n and does not 
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contribute to the derivatives. Differentiating yields 

^ = 2hM.H f n cos \n(0 - 6,)\d6 (16a) 

ar„ .'„ 



and 

B 2 E„ 



= 2hM.H [ cos [n(6 - 0„)] cos [m(0 - m )] dB (16b) 

•'o 



dr n dr 
with analogous expressions for 

dE H /Bd n , d 2 Ej,/dr n dd m , and d 2 E H /B6 n dd m . 
Evaluation of equation (16) for r b (6) = r yields 
fdE B 



, - AirrMLH, (17a) 

\dr /o 

^A = 4:tJiMM, (17b) 

d 2 E H 



= 2-KhM„H, n ^ 1, (17c) 



and all the other first and second derivatives of the applied field inter- 
action energy are zero. 

3.3 Derivatives oj the Internal M agnetostatic Energy 

The formal expression for the internal magnetostatic energy is ob- 
tained by substituting the expression for the magnetic configuration 
(6) into expression (10). In dealing with the self-interaction energy, 
it is necessary to use two coordinate systems: an unprimed system and 
a primed system. Throughout the following calculation functions of 
the spatial coordinates (r, 6, and z) are written with primes whenever 
they are of the primed coordinates. Thus M, when considered as a 
function of the primed coordinates, is written M'. The subscripted r n 
and 6„ are independent parameters and are never primed. 

The calculation begins with the evaluation of dM Jdz by differentiat- 
ing expression (6) and noting that 

£u(.r) = 5(.r) (18) 

where 8(x) is the Dirac delta function. Then 

— — = M,kg (I9a) 

oz 
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where 

k[r, r b (6)\ = 1 - 2u[r b (6) - r] (19b) 

and 

After changing the order of integration, the expression for the internal 
magnetostatic energy becomes 

e u = \m\ r r r r r r ^^^^/^v/wr'. (20) 

2. Jq Jo »o "O "—00 J— 00 S 

The factor g(z)g(z')/s contains the z and z' dependence of this integral. 
From expression (19c) it can be seen that this factor consists of four 
terms. Application of the transformation (z, z') — » ( — z, —z') to two of 
the terms under the integral sign combines these four terms into two 
terms. Making the transformation (z, z') — > (z, z), where 

z^z-z', (21) 

on the remaining terms and carrying out the integration over z yields 
the expression for the internal magnetostatic energy in terms of an 
integral over surface magnetic charges. This expression is 

»« /.cc «2t /»2x 7.7./,.,./ 

e m = m:z\ / / ^-dedd'drdr' (22) 

Jq J a •'it •mi S 

where Z is an operator defined by 

Z\ \ = f dz\h(z) - B(z - h)]\ | (23) 

s 2 = r 2 + r" - 2rr' cos (0 - 0') + z, (24) 

The factor kk' contains the r„ and n dependence of the integral so 
that the derivatives of E M may be calculated by replacing this factor 
by its derivatives under the integral. Evaluating the first derivatives 
yields 

dr„ dr h (0 7 ) dr„ + ' dr t ,{d) dr n 

= k ~ cos [n(0' - 6,,)] + fc' j^ cos [n(0 - 0„)] (25a) 

dr h or b 
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dW _ dk' dr h (e') ,, dk dr h (d) 
dd n ~ dr b (d') dd n + dr n (d) dd n 

= -k ^nr n sm [n(d f - n )] - k'-^nr n sin [n(d - 6 n )]. (25b) 
dr b or n 

Substituting these derivatives into the integral and exchanging the 
primed and unprimed r and 0. The first term becomes identical to the 
second. The derivatives of the internal magnetic interaction energy are 
then 



fr. Jo Jo Jo Jo s dr b 



dEj 

■ cos [n(0 - d„)]r'r dO' dr' dO dr (26a) 



2iji /•» r 2T r x r 2 * 



dr n dr r 



noo p2x /•«! nil -J 

= 2M\Z - 

Jo "0 "0 •'0 S 

•{*' ^ cos He - e n )] cos [m(e - e„,)] + * ^ 

•cos [n(0 - 0„)] cos W - O n )]\r'r dO' dr' d6 dr (26b) 

with analogous expressions for 

dE M /dd n , d 2 E M /dr n dd m , and d 2 E AI /dd n dd m . 

For circular domains (r b = r ) the factors k, k', dk/dr b , and dk'/dr b 
are independent of and 0' so that the integrands are periodic in 0' with 
periodicity 2tt. The range of integration of 0' may therefore be changed 
from [0, 2?r] to [0, 2tt + 0] so that after making the transformation 
(0, 0') -> (0, f) where 

r s 0' - 0, (27) 

the range of integration of both and f is again [0, 2t]. Note that now 
s 2 = r 2 + r ' 2 - 2rr' cos f + z 2 (28) 

depends only on f . 

Using trigonometric identities, the integrands of the integrals for 
the various derivatives are written as a sum of terms each of which is 
the product of a factor depending only on and a factor depending on 
f. Carrying out the integration over yields, for r b (6) = r„ , 

(*ek\ = ^m\i p r r - k> ? r>r <x (,r ' *» < 29a > 

V 3r„ /„ Jo Jo Jo s dr b 
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f w = 4^/; z r f r i (*• § + f f w « *- *, cm) 

dr /„ ./„ J ./„ 8 \ dr, ( dr b drj 

( SUM*) =1^*2 r r r- 

\ dr 2 n /„ 2 J J .'o s 

• ( fc' ^4 + ^ |- cos JiAr'r r/f dr' dr, n > 0, (29c) 

\ dr b dr b dr h I 

while all the remaining first and second derivatives are zero. Note that 
by inspection of these integrals and the definitions of k, k', and r b that 



d (BE, 



d 2 E, 



dr /a dr \ dr K 



(30a) 



and 



d*E M 

dri 



. i a. («*) _ i irMU r r r fi g a - «»»» 

2 dr \ 3r„ /„ 2 J J„ J dr b dr b s 



■r'rdtdr' dr, n > 0. (30b) 
Noting that from expressions (18) and (19b) 
dk 



Br b 



= -28(r - r„ 



(31) 



and using the definition of the Z operator given in expression (24), 
expression (29a) may be integrated with respect to r and z, and the 
second term of expression (30b) may be integrated with respect to r, ?■', 
and z. The result after some rearrangement is 



dr /„ 



= -(2Th 2 )(4TrM:)F(2r /h), 



~K~J~I = -(4xA)(4xilf ,) . ttt 

dr Q / d{2r /li) 



(32a) 
(32b) 



d 2 E„ 
drl 



= -(2^0(4^/;) 



6F(2r /h) 

d(2r /h) 



+ (h)(4rMl) "f 



LM^J)~L M 



(32 c) 



where renaming r' to r and using expression (19b) 



zr. 



' l TJ S ^ [2B(r ° ' r ° ' A) ~ 25(r ° ' r ° ' 0) 



-£(r , oo, JO +5(r , oo,0)], (33a) 
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B(r ,r f ,z)=[ f * (p 2 + zVV dr tff, 

•"0 •'0 

r 2 — 2r r cos f , 



P = >'o 



(33b) 
(33c) 



and where 

L n [(z/2r ) 2 ] = f [(z/2r ) 2 + 1(1 - cos f)] _1 (l - cosnf) df. (34) 

The L n functions are reduced to standard elliptic integral form, and 
power series expansions are obtained for both large and small values of 
the argument in Appendix A. The B function is integrated once after 
displacing the origin of the cylindrical coordinate system from to 
0' as is shown in Fig. 2. The transformation connecting the (r, f) and 
(p, ip) coordinate systems is 



p sin <p = r sin f 
p cos if = r cos f — r„ . 
After the transformation 



B(r , r, , z) = 



L h (P 2 + * 2 ) J,?0<,; 



(35a) 
(35b) 

(36a) 



f r r ' T P ^£%,r =r f . (36b) 

Jr/2 Jo {P + Z ) 

Equations (36a) and (36b) are integrated to obtain, in either case, 



B(r ,r, ,z) = [ (pi + z 2 )* d<p - \ \ z \ d<p 



(37) 



where p b is the value of p along the boimdary r = r f . For r, = ?' , p b = 
— 2r cos cp so that 




Fig. 2 — The (f, r) and (/>, <p) coordinate systems. 
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B(r , r , h) - P |/r + (2r ) 2 cos 2 *]« ^ - | I h | (38a) 
■'1/2 z 

and 

B(r , r , 0) - 2r„ . (38b) 

The remaining two terms of equation (33a) must be evaluated as a limit 

lim [B(r , r, , 0) - B(r , r f , h)] 

= lim / - [(pi + /( 2 )- - p,,] d<p + r I h I 

r^co .'0 ( 3Q) 

= 7T I /l I 

since p b approaches infinity when r, approaches infinity. Combining 
these results yields 

9 f r w/2 

' ,2 .15 



7T 

and 



F(2r Q /h) = - (2r /hf< [(h/2r f + sin 2 ^ d<p - 1 (40) 



y^ = (V2r ){2F(2roA) - - [*" [(h/2r f + sin 2 ^ A- (41) 
d(2r //i) { ir J J 

Appendix B lists the standard elliptic integral form of the force 
function F and power series expansions for large and small values of 
the argument. In Fig. 3 the force function is plotted as a function of the 
domain diameter measured in units of the plate thickness 

d/h = 2r /h, (42) 

The stability functions S„ , also shown on this plot, are defined in Sec- 
tion 5.1. 

IV. THE ENERGY VARIATION — ORIGIN OF TERMS 

Summing the results of the last section according to expression (7), 
the total energy variation expression (11) is 

AE = [27r/i<r„, + ±wr hM,H - (2irh 2 )(4ir]\l])F(2r /h)] Ar a 
+ \ [iMIM - (4x/0(47r.l/!) f|^](Aro) 2 

+ 1 1 fe *"' + «"■* - <**<«* §8? 

+ (h)(^M]) I 2 [L„((V2r ) a ) - L„(0)]}(Ar n ) 2 + 3 (43) 
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Fig. 3 — The magnetostatic radial force function F and stability functions. 
So — Sio, Si = 0, as functions of domain diameter to thickness ratio, d/h. 

where F is defined by expression (33) and plotted in Fig. 3, the L n are 
defined by expression (34), and all terms not explicitly stated are equal 
to zero. The remainder of this section treats the physical origin of the 
terms in the energy variation expression (43). 

4.1 The Generalized Forces 

The coefficients of the linear variation terms are the negatives of the 
generalized forces. All forces except the r force are identically zero, 
which for a circular domain is a consequence of the rotational symmetry 
of the system. The first term in the coefficient of Ar„ is the product of 
the wall energy density <r w and the rate of change of wall area with re- 
spect to r , 2irh. The second term is the product of the external field 
interaction energy density 2M,H and the rate of change of domain 
volume with respect to r , 2whr . The third term is the rate of change 
of the internal magnetostatic energy with respect to r . 

The internal magnetostatic force may be identified in expression (43) 
and using expressions (32), (33), and (39) may be written in the form 



_(dE JL ) 
\dr J 



= 2irh 2 (±TrM 2 .)F(2r /h) 



= (2vroh)(2M.)\4irM. 



+ 



4M. 



[ r r-L±M - r r 

Jo h (p"+W Jo Jo 



r dr d£ 



P . 



(44a) 



(44b) 
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where p 2 is given in expression (33c). In expression (44b) the first factor 
in parentheses is the domain wall area, the second is the change in 
magnetization at the wall as the wall moves, and the quantity in braces 
has the form of an H field whose origin will now be interpreted by super- 
position of sources. The internally produced field arises from the super- 
position of the internal field of a plate uniformly magnetized normal 
to its surface with magnetization magnitude M, and two disks of mag- 
netic charge of uniform magnetic surface charge density ±2M, and 
radius r . The first term within the braces is thus the demagnetizing field 
of the infinite plate of uniform magnetization. The second term is the 
difference in magnetostatic potential between a point on the edge of a 
disk of magnetic surface charge of uniform density 4M, and a point 
removed a distance h from this point in a direction normal to the plane 
of the charge disk divided by the distance h. This is just the z-averaged 
z-component of the field produced along the wall by the two charge 
disks since 

where ft denotes the magnetostatic scalar potential and z is measured 
from the edge of the disk. Comparing expressions (44a) and (44b), the 
total internally produced 2-averaged z-component of the magnetic field 
along the domain wall is 

<ffir.>- = -(±*M.Kh/2r )F(2r /h) (46) 

so that the total force per unit wall area (averaged over z) is 



-Mf;)r-°t- 2 ^ H + ^- 



(47) 



The first term is the product of the wall energy density and the wall 
curvature and always corresponds to an inward directed force. The 
second term is the change of magnetization at the moving domain wall 
times the z-averaged z-component of the total field at the wall. [The 
problem may initially be set up using this fact (Ref. 2, pp. 1922-1925).] 
The properties of the force function will now be examined in some detail. 
From expressions (43) or (44a) the first order variation in internal 
magnetostatic energy, when ?- is varied, is 

AE U = -2(T}i*K±irM 2 .)F(2r /h)Ar . (48) 

The plot of F in Fig. 3, and the expansions for large and small values 
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of the argument show that the force function is everywhere positive, 
is monotonic increasing, and has a negative second derivative. Since 
the force function is everywhere positive, the internal magnetic inter- 
action energy at all times acts in such a way as to expand the domain. 
Section 4.2 treats the effect of the slope and curvature properties of the 
force curve on domain size and stability. Substituting the expansion of 
the force function for small values of the argument (138d) into expres- 
sion (48) produces the energy variation for small values for r /h, 

AE„ = {2irhr 2MX-4*M a ) + (47r-M 2 )16r 2 

- (7r/i 2 )(4rM 2 .)2[K2roA) 3 - A(2r A) 8 + ■ • ■]} Ar . (49) 

(In the remainder of this section frequent reference will be made to the 
properties of F and the L n given in Appendices A and B.) The inter- 
action of the magnetization with the existing field from the infinite 
dipole sheet produces the first term in expression (49). This may be seen 
by comparison with expression (47) and observing that the field in- 
ternally generated in the infinite dipole sheet with no reversals is 
— 4rAf, . In Fig. 3 a dashed line through the origin with numerical 
slope one represents this term and forms the small r /h asymptotic of F. 

The second term in expression (49) is the only thickness independent 
term in the expansion and therefore must be identical to the variation 
of self-energy of the two disks of magnetic charge which form the ends 
of the reversal when r is varied. Since the interaction with the infinite 
charge sheet and the self-energy of the disk have been taken into ac- 
count, the remaining terms are the mutual interaction of the magnetic 
charge disks. 

For large r /h, an energy expansion in terms of h/r is appropriate. 
Substituting the expansion of the force function for large values of the 
argument (138c) into expression (48) yields 



AE M = -/i 2 (47ril/ 2 ){[l + Mh/2r f + 0,] 



+ [2 - \{h/2r ) 2 + 4 ] ln 



«* 



Ar ■ (50) 



This expansion obscures the identity of both the infinite sheet magnetic 
field term and the charge plate self-energy term so that a local (to the 
wall) magnetic energy lowering per unit line length description appears 
appropriate. However, the energy reduction per unit line length to 
lowest order in 2r Q /h is 
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E ^f (domain) — E M (uniform magnetization) 



- (4tM 2 ,) In 



42n, 
e* h 



(51) 



so that the energy lowering per unit line length for the domain of in- 
finite diameter is infinite. [Equation (51) is obtained by integrating 
equation (50) to lowest order. The integration constant is determined 
to be zero by term by term integration of the expansions of F for large 
and small values of the argument and comparing at 2r /h = 1.] The 
conclusion that the energy lowering per unit line length for an isolated 
straight line reversal may also be obtained by considering the energy 
lowering in a strip reversal when the strip width approaches infinity. 
The author's intention at the outset of this entire calculation was to 
calculate the numerical value of this magnetic energy reduction per 
unit wall length. The internal magnetic interaction, however, retains 
just enough of its global character when the domain is very large so that 
no finite limiting value for this energy reduction exists. 

The internally generated magnetic field at the wall of the domain, 
for large r /h is obtained from expression (46). To lowest order it is 






h 



(52) 



which approaches zero as the diameter approaches infinity as it must, 
since for an infinite straight line magnetization reversal, symmetry 
requires that the z-component of the field be zero along the reversal. 

4.2 The Stiffness Matrix 

The second variation of the energy with respect to the Fourier coeffi- 
cients describing the domain determines the stability of the domain. 
Since the stiffness of the domain with respect to externally applied 
forces is proportional to the coefficient of the bilinear form which is the 
second variation of the energy, the matrix formed by these coefficients 
is called the stiffness matrix. The stiffness matrix is composed of three 
independent submatrices. The second derivatives of the energy with 
respect to the Fourier amplitudes form the radial stiffness matrix; 
the second derivatives of the energy with respect to the Fourier phases 
form the angular stiffness matrix; and the derivatives of the energy with 
respect to one Fourier amplitude and one Fourier phase form the mixed 
stiffness matrix. The derivative of the energy with respect to r„ and r m 
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are called the (n, in) radial stiffness matrix element, with similar nota- 
tion for the other submatrices. 

All derivatives not explicitly exhibited in expression (43) are zero. 
Thus, the angular stiffness matrix and the mixed stiffness matrix are 
zero and the radial stiffness matrix is diagonal so that the system is 
completely metastable with respect to angle and the amplitudes are 
normal modes of the system for small amplitudes. 

The (0, 0) radial stiffness matrix element is simply the derivative of 
the negative of the radial generalized force so that no further discussion 
of it is necessary. It should be noted that the derivative of the internal 
magnetostatic term with respect to wall position is not directly related 
to the radial field or potential at the wall since the derivative used in 
computing the radial field at the wall must be taken with the wall posi- 
tion held fixed. 

4.2.1 The Radial Stiffness Matrix Elements for n ^ 1 

The diagonal radial stiffness matrix elements, for n ^ 1, are the sum 
of four terms in expression (43). The first term, which always has a 
stabilizing effect, is the increase in total wall energy due to the lengthen- 
ing of the wall caused by the deviation from a strictly circular shape. 
Imposing a sinusoidal variation of amplitude As onto a straight line 
produces a relative increase in length of 

S \ A n / 

The corresponding wavelength in expression (43) is 

X„ = 2=&. (54) 

n 

The wall energy term in expression (43), 

AE Wn - v,airr h(^=^y , (55) 

is thus the product of the wall energy density, the wall area, and the 
variation in wall area per unit area. Notice that the relative variation 
in wall length or area is independent of the wall curvature, l/r , to 
lowest order in the amplitude of the variation. 

The second order change in volume of the domain interacting with 
the externally applied field produces the second term in the radial stiff- 
ness matrix elements while the rate of change of the internal mag- 
netostatic forces at the wall produces the third term. The sum of the 
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second and third terms is one-half the (0, 0) radial stiffness matrix 
element. This factor of one-half relates to the fact that a variation of 
Ar n , n ^ 1 produces only one-half the mean square variation r b (6) 
as is produced by an equal variation in r . This shape-independent, 
second-order variation in energy arises from the variation in the gen- 
eralized forces, or fields at the wall, when the domain radius is varied. 
[See also the steps leading to expression (30b).] 

4.2.2 Translation Invariance 

The requirement of translation invariance in the infin ite plate com- 
pletely determines the (1, 1) radial stiffness matrix element. Consider 
a cylindrical domain of radius r with a cylindrical coordinate system 
placed at its center. Under a displacement of the coordinate system of 
magnitude s in the 8 = -n- direction, the description of the boundary 
in the new coordinate system is 

1 2 1 <r 

r h (0) = ro - - i- + s cos d + - - cos 20 + 0., . (56a) 

4 r 4 r 

Thus, to second order in s, term by term comparison with definition 
(3) yields 

Ar„ ■ — t— , Ar, = s, and An = j— ■ (56b, c,d) 
4 r t r 

The formal change in energy under this displacement (11) is 

**- ®.(-3 +!©/+*■ 

Obtaining (dE/dr ) and (d 2 E/dr l ) from expression (43), and substitut- 
ing expressions (84), (85), (86), (100), and (138a) verifies that 

d'E T \ 1 (dE T \ /(-n\ 



dr\ n 2r \dr 

The coefficient of s 2 in expression (57) is thus zero as required by transla- 
tion invariance, and further the (1,1) stiffness matrix element is zero 
whenever the total radial generalized force is zero. 

4.2.3 The Magnetostatic Stiffness Terms 

The interpretation of the radial stiffness matrix elements for the 
higher n values is now considered. As in the case of the generalized 
forces, examination of the expansions for small r allows the self-inter- 
action energy of the two charge disks which make up the ends of the 
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domain to be separated from the mutual interaction of these charges. 
The variation in the internal magnetostatic energy due to a variation 
in some r n for a circular domain is in general from expression (43) 



AE Un m 47ril/: 



-^^f^WS-HW 



ft £ 1. (.59) 



Separating the h independent and h dependent terms of the power 
series representation in powers of 2r /h uniquely separates the above 
expression into two parts, one part representing the self-interaction of 
the charge disks and the other representing the mutual interaction of 
these disks. The h independent terms then represent the self-interaction 
forces of the charge disks and the h dependent terms represent the 
mutual interaction forces. In the expansion of L„ for large (h/2r ) 2 , 
expression (129a), all terms of L n [(h/2r ) 2 ] are h dependent. Using the 
large (h/2r ) 2 expansion of F, expression (138d), and the expressions 
for Z/„(0), (115) and (116), the thickness independent part of expression 
(59) is 

A#A,„(Self) = 4*MVro(* - 4 E ^pr^Ar,,) 2 , n 2 1. (60) 

This energy variation contains a term which results from the variation 
in the overall size of the disks of charge as well as the shape dependent 
terms. The size variation term will now be identified and subtracted 
out so that the shape dependent part of the self-interaction energy may 
be seen explicitly. From expression (49) and the discussion following it, 
the ratio of the variation in energy of two isolated disks to the variation 
in disk area is (47rilf*)(16r /27r). The variation in disk area for a varia- 
tion in r n for n ^ 1 is (ir/2)(Ar„) 2 so that the change in self -energy of 
the two disks, other than that due to their mutual interaction or change 
in overall size, is 

AE iIn (Self-Shape) 

0, n = l (61a) 

-(4xil/ 2 )4r (|:^ J n -)(Ar n ) 2 , n > 1. (61b) 

It is not surprising that the variation of ?'i produces no shape related 
energy change, since from expression (56) this variation is to lowest 
order a displacement with a size change coming in second order. It is 
seen that the terms which remain after cancellation all come from L„(0). 
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In expression (59) the first term is independent of n and the 
— 47rM^ , L„(0)(A? , n ) 2 term has been identified with the variation in the 
self-energy of the charge disks. The term 4:vM 2 .roL n [(h/2r ) 2 ](Ar n ) 2 must 
therefore contain all of the shape dependent part of the charge disk 
mutual interaction energy. This term also contains a contribution due to 
the variation in the total amount of charge and contribution due to the 
shape independent, general smearing out of the charge distribution. 
Since the second order change in the total amount of charge is inde- 
pendent of n for n ^ 1, these two contributions may be removed from 
the mutual interaction energy variation by replacing L n by L n — L x . 
The remaining mutual interaction energy variation is specifically due to 
the shape of the variation. This energy variation is 

A#A,„(Mutual-Shape) 



h 



(Ar n ) 2 , n ^ 1 



= (47tM> [l b (^) - U 

n ^ 1 (62) 

where the final form is obtained using the expansion for L„ , equation 
(131), and the M n _ m are the constants of the expansion. The interaction 
energy of planar multipoles of order n and higher has the form of equa- 
tion (62), as it must since the variation in the charge distribution for 
each n may be expressed in terms of such multipoles. 

The variation in internal magnetostatic energy due to a variation of 
r n , in the infinite sheet, for large 2r /h, to lowest order in h/2r , is 

AE Mn = (4xil0(^) 

• [-2n 2 In | 4 ^ | - 2n 2 - 2 + (4n 2 - 1) £ ^pn](Ar„) 2 , 

n ^ 1 (63) 

using equation (59), the large 2r /h expansion of F (138c) and of 
L„(/i74r 2 ) - L n (0), (105), (116) and (125). 

The charge-disk self-interaction energy is not evident in this ex- 
pansion because it is exactly cancelled by the leading term of the mutual 
interaction energy. In contrast to the energy reduction per unit line 
length for a straight line reversal in an infinite sheet (which has no 
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finite value), it is possible in the case of this variation of the domain 
structure to compute the energy variation per unit line length. In terms 
of the wavelength of the variation X, defined in expressions (54) or 
(127) and in the limit of r — » °° , the total variation in energy when r„ 
is varied is [using expressions (43) and (55) and the limit (128)] 



Aff, 

27JT,, 



j[xV,„//i - (4rJfftr In | ^ ] £ + 0,(^)}(AO 2 . (64) 

Comparison of expression (64) with expression (53) shows that the 
magnetostatic energy variation per unit line length for a circle of infinite 
diameter is the product of the magnetostatic energy density constant, 
the variation in line length, and the logarithm of a maximum effec- 
tive interaction distance, 4eX/x. (The maximum effective interaction 
distance for the magnetostatic energy lowering per unit line length is 
proportional to r .) Hagedorn has computed the magnetostatic energy 
variation per unit line length for the case of a sinusoidal variation im- 
posed on an infinite straight line reversal. The calculation was carried 
out by considering the energy variation produced by a sinusoidal applied 
to a strip domain pattern in the limit of infinite strip width. The result 
of this calculation is 

LEul (unit length) = -(4rM*)rln | X/(2.111/i) | (/i/X) 2 (Ar) 2 , (65) 

which differs from the result for the infinite circle by the constant inside 
the logarithm. 

4.3 Summary 

The physical origin of terms of the energy variation has thus been 
traced in the limiting cases of both large and small r /h. In either of 
these limiting cases, it is thus possible to develop intuition with regard 
to the behavior of the domains. Since, as has been shown, the inter- 
pretation of the meaning of the energy terms in the limiting cases is 
qualitatively different, the development of intuition in the transition 
region is quite difficult. In many device applications this transition 
region is the preferred region of operation, making the use of analytical 
and numerical methods a necessity. 

V. THE SIZE AND STABILITY OF CYLINDRICAL DOMAINS 

The energy variation expansion (43) in principle contains all cylin- 
drical domain size and stability information. This section treats briefly 
the use of this expression in the determination of domain size and 
stability. The only non-zero generalized force in expression (43) is the 
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uniform radial force. When this force is set equal to zero (the force 
equation), the system is in equilibrium. Thus the condition that the 
system be in equilibrium provides, given a material and plate thick- 
ness, an equation relating domain size and the applied field. The loca- 
tion of the zeros in expression (43) (all terms not explicitly exhibited 
are zero) shows that the system is completely metastable with respect 
to angle and that the radial stiffness matrix is diagonal. The radial 
amplitudes are thus quasi-normal modes, and the study of stability 
reduces to the study of the stability of the individual radial amplitudes. 

5.1 Normal Form o] the Energy Expansion 

Before proceeding with the discussion, it is appropriate to introduce 
some new notation and to rearrange the energy variation expansion 
into what will be called normal form. Since the stiffness matrix is of 
interest only when the domain is in equilibrium, the applied field H 
is eliminated from it using the force equation. The geometrical de- 
pendences of the various magnetostatic stability terms are then com- 
bined and normalized to the wall stiffness term by defining the "stability 
functions" as 

S (d/h) - F(d/h) - dj- d F(d/h) (66a) 

and 
S n (d/h) - -^ZTl {^(d/h) + ^ (d 2 /h 2 )[L n (h 2 /d 2 ) - L n (0)]}, 

n ^ 2. (66b) 

The Si function is undefined or may be taken to be zero since transla- 
tion invariance in the infinite plate requires that the (1, 1) stiffness 
matrix element the identically zero whenever the generalized radial 
force is zero, as is assumed to be the case here. The S n functions are 
plotted in Fig. 3 up to S 10 ; they are given in standard elliptic integral 
form together with power series expansions for large and small values 
of the argument in Appendix B. The domain diameter, d = 2r repre- 
sents domain size in this section. The normal form of the energy expan- 
sion is written as a function of the ratios of the three fundamental lengths 
of the system: the plate thickness h, the domain diameter d, and the 
"characteristic length" defined by 

' - i^r <67) 
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The characteristic length depends only on the type of material used. 
Dividing the energy variation expansion (43) by the normalizing 
energy 2{±TrM])(irti i ) and introducing the notation of the preceding 
paragraph, the normal form of the energy expansion results: 



A-E, 



2(47n'l/ L ;)(7r/i 3 ) 



I d H 

h + h 4ttM, 



-4 



h-1-H 



i - <i 



Ar„ 



+ £ (»' - i) 



< 



Ar B 
h 



+ 3 



(68) 



In expression (68) the coefficient- [l/h+ (d/h)(H/4*M t ) - F(d/h)] 
is the normalized radial force. Setting this force equal to zero yields 
the normalized force equation. The remaining bracketed quantities 
[l/h — S n (d/h)] are proportional to the diagonal elements of the stiff- 
ness matrix, and are called "stability coefficients." For uniform radial 
variation, the stability coefficient has the opposite sign from the (0, 0) 
element of the radial stiffness matrix; thus this stability coefficient 
is negative whenever the domain is stable. For the other r n variations, 
on the other hand, the stability coefficient has the same sign as the 
corresponding element in the stiffness matrix, and these stability co- 
efficients are positive whenever the domain is stable. 

5.2 Graphical Solution of the Force Equation 

A graphical solution to the force equation 



1 4- i -JL- - w(«\ = 
k^~h4irM. \hJ 



(69) 



may be obtained by constructing a straight line on Fig. 3 whose inter- 
cept with the vertical axis is l/h and whose numerical slope is H/4irM, . 
The intersections of this straight line with the F curve are then the solu- 
tions to the force equation. 

As was stated in Section 5.1, (i) the force function has a positive first 
derivative and negative second derivative for all nonzero values of its 
argument, (ti) it is zero and has a first derivative of unity when its 
argument is zero, and (Hi) it becomes logarithmic for large values of 
its argument. From these properties and examination of Fig. 3, several 
properties of the solutions to the force equation may be appreciated. 
For negative values of the applied fields, there is only one solution to 
the force equation. Examination of the sign of the radial force which 
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results when the diameter is varied about the solution diameter while 
all other variables are held fixed shows that this solution is unstable. 
For small positive applied fields, there are two solutions to the force 
equation, the larger diameter solution being radially stable, the other 
radially unstable. However, a radially stable solution does not guaran- 
tee that the system is stable with respect to all possible deformations, 
and this must be investigated separately. As the applied field is increased, 
the two solutions move closer together until they coalesce. When the 
applied field is increased beyond this point, there are no solutions. 
Since the function F is asymptotic to a straight line through the origin 
having unit slope, the solutions will always vanish for a value of the 
applied field which is greater than 4irM. . Stable isolated cylindrical do- 
mains thus exist only in the presence of an applied field having magni- 
tude between zero and 4:irM, and polarity tending to collapse the 
domain. 

5.3 Graphical Determination 0/ Domain Stability 

The stability coefficients are determined graphically by constructing 
a horizontal fine at height l/h on the force stability graph. Metasta- 
bility for each normal mode of deformation occurs at the intersection 
of this line with the corresponding stability function. Since the stability 
functions are monotonic, the diameter of metastability of each normal 
mode of deformation is uniquely defined and forms the boundary be- 
tween the regions of stability and unstability. The circular domain will 
be stable with respect to all variations when its diameter is greater than 
the radial metastability diameter and less than the metastability diam- 
eter for a variation with a rotational periodicity of two. The normal 
variations with rotational periodicity two are referred to as "elliptical" 
deformations. When the domain is stable with respect to elliptical 
deformation, it is necessarily stable with respect to the variations 
of higher spatial frequency since the stability functions of higher spatial 
frequency lie progressively (with respect to n) below the elliptical sta- 
bility function. The radial stability function S and the elliptical stability 
function S 2 thus form the boundary of the region of total cylindrical 
domain stability. Therefore, given the magnetic material type and 
plate thickness, the range of stable domain diameters and the correspond- 
ing applied fields may be determined with the aid of these functions. 

5.4 Minimum Domain Diameter 

For any given value of l/h, the minimum domain diameter is the 
collapse diameter determined by S . The domain diameter measured in 
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units of the characteristic length is d/l = (d/h)/(l/h) which is the in- 
verse of the numerical slope of a line drawn in Fig. 3 from the origin 
to the operating point. The line of maximum slope, which both passes 
through the origin and contacts the *S curve at at least one point, thus 
determines the smallest domain diameter attainable in a given material. 
The coordinates of this contact point are d/h ^ 1.2 and l/h « 0.3, so 
that the minimum attainable domain diameter is 

d mia » U. (70) 

VI. RANGE OF VALIDITY OF THE MODEL AND THE QUALITY FACTOR 

At the present time, no quantitative evaluation of the range of valid- 
ity of the domain structure model used here has been carried out. The 
qualitative discussion given here, it is hoped, will provide the reader 
with an appreciation of the magnitude of the effects produced by the 
relaxation of the various constraints artificially imposed by the model 
and the dependence of these effects on the system parameters. It has 
been assumed that domain walls are cylindrical, have zero width, and 
have a definite energy per unit area which is independent of wall orienta- 
tion or curvature, and that the magnetization lies perpendicular to the 
surface of the plate. Section 6.1 treats the effect of the relaxation of 
the cylindrical wall approximation only. In Section 6.2, the other 
assumptions are all shown to be coupled using the simplest uniaxial 
material model. A single dimensionless material parameter q, which 
complements the characteristic length I in characterizing circular domain 
materials, is used to express the results obtained from the simplest 
material model. 

6.1 The Cylindrical Wall Approximation 

The discussion of the cylindrical wall approximation uses the coordi- 
nate system and domain configuration of Fig. 1 except that the walls 
are allowed to curve as shown in Fig. 4. The radius function, r b {d, z), is 
determined by the requirement that it minimize the total energy. The 
Euler equation which results from this two dimensional field variational 
problem is an integro-differential equation similar to those which appear 
in Hartree self -consistent field calculations. No solution of this equation, 
numerical or otherwise, has been attempted or is contemplated at the 
present time. The Euler equation consists of terms arising from: the 
wall energy, the interaction of the magnetization with the applied field, 
the self-interaction of the magnetostatic charges at the surface of the 
plate, the self-interaction of the charges produced by the slope of the 
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Fig. 4 — Cross section of a noncylindrical, near-circular, domain. 



domain wall, and the mutual interaction of the plate surface charges 
with the domain wall magnetic charges. Boundary conditions (obtained 
from the appropriate transversality condition 7 ) require that the wall 
surface be perpendicular to the plate surface at all intersection points 
(3 in Fig. 4). Physically (since in the model used here the crystal is 
assumed to be strain free) the surface cannot interact with the domain 
wall, and therefore the wall must intersect the surface at right angles. 

Although it is not clear that domains having a roughly conical shape 
are ruled out, it will be assumed that the domain has reflection symmetry 
through the central plane of the plate and that the radius is a function 
of z only, r b (z). In this case the wall must be vertical at the central 
plane as indicated at Q in the figure so that the single parameter b 
represents the magnitude of the wall bulging. Since the Euler equation 
requires the curve to be smooth, there must be an inflection point, S, 
between 6 and 3. The wall area, and thus total wall energy, is a quadratic 
increasing function of the wall curvature so that the concentration of the 
curvature at the center and ends of the wall, produced by the trans- 
versality and symmetry conditions, tends to reduce the wall bulging. 

The radial field at the domain wall from the charges at the surface 
of the plate is directed as shown in Fig. 4. The effect of the interaction 
of the magnetostatic charges due to the slope of the wall with the radial 
component of the field from the surface charges is destabilizing for 
either positive or negative bulging. This interaction produces a negative 
quadratic term in the total energy. However, at the plate surface, where 
the magnitude of the radial field is greatest, the transversality condition 
requires that the charge density produced by the wall slope is zero so 
that the magnitude of this negative term is small. The z component of 
the field from the charges on the surface of the plate determines the 
direction of bulging. (The applied field, being uniform by assumption, 
need not be considered.) Along an initially cylindrical wall the internal 
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field is everywhere directed, so as to make the domain expand, and 
attains its greatest magnitude at the center plane of the plate. It there- 
fore provides a linear term in the total energy which tends to bulge the 
wall in the positive direction as shown in Fig. 4. 

Thus, for near cylindrical walls, the bulging is determined by the 
interaction of this force (tending to bulge the wall) with the wall energy 
(acting to stabilize the wall) and the radial field (acting to destabilize 
the wall). The self -interaction of the wall charges enters only as a 
higher-order term. It should be noted that the transversality condition 
acts both to strengthen the stabilizing term and weaken the destabilizing 
term. 

The relevant dimensionless wall energy for the wall bulging problem 
is l/h = (T w /(h4:TrM 2 B ). Wall bulging is expected to decrease with increas- 
ing wall energy. A second independent effect related to l/h may be 
appreciated by inspection of the S and S 2 curves in Fig. 3. It can be 
seen from Fig. 3, equation (68), and the discussion following it that, 
since the S and S 2 curves bound the region in which stable circular 
domains exist, d/h must increase with increasing l/h. By symmetry, 
the ^-component of the internally generated magnetic field at a cylindri- 
cal wall is zero for a domain of infinite diameter and clearly increases 
monotonically as the domain diameter to thickness ratio decreases. 
Thus, as the plate is made thicker, the bulging force becomes stronger 
and the stabilizing force becomes weaker. Since several independent 
effects cooperate to increase bulging with increasing plate thickness, 
the onset may be quite rapid when it does occur. Domain collapse data 
taken at d/h ~ 1 is in good agreement with predictions made on the basis 
of equation (68) and Fig. 3. 8 This then provides some indication that 
the cylindrical wall approximation remains valid at this thickness. 

6.2 The Quality Factor 

The discussion of the approximations other than the cylindrical wall 
approximation uses a polar (ilf„, tj, v) coordinate system where rj is the 
polar angle and v is the aximuthal angle to specify the orientation of 
M , (See Fig. 5). The polar axis is taken to be the 2-axis of the preceding 
sections. The domain wall is taken to be planar with its position and 
orientation specified by a plane at its center. The axis through the origin 
in the direction of the wall normal is denoted by £. The position of the 
central wall plane is denoted by £ . The orientation angles of the wall 
normal are denoted by v w and r} w , (see Fig. 5). 

In the simplest uniaxial material whose easy axis is the 2-axis the 
magnetic energy density for a planar wall is (Ref. 9, pp. 189-192) 
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Fig. 5 — Coordinate system for specification of domain walls. 
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(71) 

where A is the exchange energy density coefficient, K u is the anisotropy 
energy density coefficient, H is the externally applied field, and the last 
term is obtained by integrating V • B = 0. For a uniformly magnetized 
material in the absence of applied or internal fields, this expression 
reduces to p E = K u sin 2 -q which has absolute minima at rj = and 
v = 7r. The 2-axis is thus the easy axis as is required for consistency 
with the preceding sections. 

The anisotropy energy density coefficient is sometimes expressed in 
terms of the effective anisotropy field H a m 2k u /M,. The quality factor 
is now defined as the dimensionless anisotropy energy coefficient or 
dimensionless anisotropy field 



Q = 



K„ H a 

2wM': -iwM, 



(72) 
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6.2.1 The Nucleation Field 

When a bias field H is applied in the positive z-direction and demag- 
netizing fields are neglected the energy density is 



PE 



= K„ sin 2 77 — HM, cos r? 



(73) 



which, for H < H a , has a local minimum at magnetization orientation 
tj = 7T and an absolute minimum for 77 = 0. When H > H a , only the 
minimum at tj = remains. In a perfect crystal the effective anisotropy 
field is thus the field at which the magnetization becomes unstable 
with respect to reorientation (assuming it is initially oriented in the 
negative z-direction). If a reorienting field is applied locally (local but 
over a region whose dimensions are much greater than a wall width so 
that the effect of exchange forces can be neglected), then H a is the total 
local field required for the nucleation of a domain at that locality. If 
the nucleation field H N is understood in this sense, then in a perfect 
crystal q is the nucleation field measured in units of 4iril/.: 



H, 



4*M. 



= Q- 



(74) 



In an imperfect crystal H N /4irM, may be either larger or smaller than q. 
If it is larger, the material may be expected to have a high wall motion 
coercivity. 

6.2.2 Susceptibility 

When a transverse bias field H t {i\ = w/2) is applied and demagnetizing 
fields are neglected, the energy density becomes 



Pb 



= K u sin 2 r] — H,M, sin 7? 



which has stable magnetization orientations 



sin "' W - sin "' Gar: \i ' 



v = 



jr 

2 ' 



H, < H a ; 
H, ^ H a . 



(75) 



(7(5) 



The transverse susceptibility is therefore 
dM, j!- 

(H t * H a ) 



(H, < H a ) 



(77) 
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where M , = M, sin 77 is the component of the magnetization in the 
direction of H t . Thus, the susceptibility to tipping of the magnetization 
by a transverse field is inversely proportional to g. 

6.2.3 Wall Energy and Wall Width 

Consider now a planar Bloch wall, J[-M = 0, between two regions 
whose magnetization at points far from the wall lies along the two easy 
directions, 77 = and jt, and again assume that there are no applied 
fields or fields produced by boundary surfaces. Under these conditions, 
the magnetic configuration is determined by minimization of the wall 
energy per unit surface area which in this case is (Ref. 9, pp. 189-192) 



-£ MS" +*•*'• 



a* 

Carrying out the minimization results in 

a w = 4(AKrf 
for 

£ — &> = ~ h 1°S tan 

7T 



«• 



where 






(78) 



(79) 



(80) 



(81) 



is the wall width. The definition of wall width is somewhat arbitrary 
since the wall extends over all space. In this case, following page 191 of 
Ref. 9, it is chosen so that the magnetization would complete its entire 
rotation of t radians in a length l w if the entire rotation took place at 
its maximum rate, the rate at the center of the wall. 

The ratio of the characteristic length, equation (67), to the wall 
width is 



L 



2 

= - g, 

IT 



(82a) 



so that the ratio of the minimum domain diameter, equation (70), to 
the wall width is 



d min _ 4/ _ 8 

In L T 



(82b) 



The approximation of zero wall width thus improves as q becomes 
larger. 
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The approximation that the wall energy is independent of wall curva- 
ture is clearly related to the wall width. At large distances from the 
planar wall equation (80) becomes 

| t, - tjo | = 2exp(-7r I £ - & I /U (83) 

where i) is the appropriate equilibrium orientation of the magnetization 
at a distance far removed from the wall. Such an exponential relation 
will hold for the approach to any stable equilibrium orientation in the 
presence of isotropic exchange. The change in energy of the wall due to 
overlapping of the tails of the wall as the wall is curved is clearly related 
to q, becoming larger as q becomes smaller. In order to solve for the 
dependence of the wall energy on curvature it is necessary to solve the 
entire (including magnetostatics) micromagnetics problems. 

6.2.4 Summary 

The preceding results may be summarized by noting that the higher 
the q value, the more closely the simple uniaxial model obeys the con- 
straints of the domain model used in the previous sections. It is clear 
that, for domains of the type considered to exist at all, q must be greater 
than one. For device operation, q should probably have a value greater 
than two. 

VII. CONCLUSIONS 

The theory of cylindrical magnetic domains yields conditions winch 
predict the size and stability of these domains and provides an estimate 
of the range of applicability of the model used. The results of theory 
appear to be accurate in a range useful in the construction of circular 
domain devices. 

The domains considered are isolated right circular cylinders in plates 
of uniaxial magnetic material of uniform thickness cut so that the plate 
normal is parallel to the easy axis. The first and second order energy 
variations which result from a general small deviation from the strictly 
circular shape determine domain size and stability. The energy method 
was chosen in preference to the magnetostatic field method because of the 
uniformity it provides in accounting for the forces in both the equi- 
librium and stability problems. The integrals arising from the energy 
method are interpreted physically in terms of fields and interacting 
charges. The physical interpretation of the integrals is quite different 
in the limiting cases of very large or very small domains. The integrals 
are related to special cases of the fields of uniformly charged disks 
computed by C. Snow and tabulated by N. B. Alexander and A. C. 
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Downing. 11,12 The present work obtains the needed properties of the 
integrals (expansions, recursion relations, and others) directly from the 
definitions. 

When the energy variation is described in terms of a Fourier decom- 
position of the domain radius function, only the generalized force cor- 
responding to a change in domain size is non-zero and the stiffness 
matrix is completely metastable with respect to angle (phase) and 
diagonal with respect to the Fourier amplitudes. Since the Fourier 
amplitude stiffness matrix elements are all found to be distinct, the 
description is unique and may be described as a quasi-normal mode 
description. 

The normal mode description is summarized by a single graph from 
which many domain properties may be determined by construction. 
Cylindrical domains exist only in the presence of a bias field directed 
so as to tend to collapse the domains and having a magnitude between 
and 4x71/, . The uniform radial collapse of the domain and the run- 
out of the domain into an initially elliptical shape bound in the region 
of stability. The minimum attainable domain diameter in a given 
material is d min ~ 4Z occurring a plate thickness of ~ 4Z. It is estimated 
that the cylindrical wall approximation begins to become doubtful at a 
plate thickness greater than 41. In order for cylindrical domains to 
exist, H a « 471-il/, and in general approximations such as the approxi- 
mation of zero wall width become more accurate for H a » 4irM . (d mia /l w = 
8HJ4:ir 2 M. where l u is the wall width). 

It is interesting to note that since stable cylindrical domains of a 
definite size exist in the total absence of wall motion coercivity and may 
be freely moved, they form a relative, easily observable, classical model 
for illustrating several particle-field concepts. They may be considered 
a two-dimensional particle which is produced as a singularity of finite 
extent in an underlying three-dimensional field (the magnetization). 
Cylindrical domains are particularly useful for demonstrating the con- 
cept of identical particles since, while it is possible to put identifying 
marks on domain locations, it is not possible to mark individual do- 
mains. (Cylindrical domains do exist in two species which may be 
distinguished by the direction of rotation of the spins in the domain 
wall. 13 All attempts to observe this difference up to the present time 
have been unsuccessful.) 
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APPENDIX A 

Integrals of Cylindrical Domain Theory 

This appendix contains the reduction to standard form of the elliptic 
integrals which arise in the theory of cylindrical domains in plates of 
infinite extent and power series expansions of these integrals. All the 
properties of the functions obtained here are used in either the physical 
interpretation of the energy variation expansion or in generating the 
numerical values of the force and stability functions. 

It is convenient to define functions U and V which appear repeatedly 
in cylindrical domain theory. The elliptic integrals which appear in the 
final results of the theory appear only in the forms U and V, U being a 
function of only the complete elliptic integral of the second kind and V 
being only a function of the complete elliptic integral of the first kind. 
Because of the form of the U and V functions, it has proven easier to 
obtain the needed properties, (such as the series expansions) directly 
from the integral definitions rather than deducing them from the tabu- 
lated properties of elliptic integrals. 

The latter half of this appendix treats the properties of the L„ func- 
tions. A recursion relation is obtained and used to reduce the L n to 
functions of U and V. Power series expansions of the L n are obtained 
directly from the definition (34) . 

A.i Definition of the U and V Functions 

The functions are defined in the alternate forms 

U(x) m [' [x + Kl - cos a)]* da (84a) 

= 2 T (x + sin 2 /3)*d/3 (84b) 

Jn 



CYLINDRICAL MAGNETIC DOMAINS 3321 

= 2 / (x + 1 - sin 2 y) h dy (84c) 

= 2(* + D^(r^) 



(84d) 



aud 



V( x ) = f [x + Ml - cos a)]" 1 da (85a) 

= 2 ( (.c + sin 2 /3)~* dp (85b) 

•'0 

= 2 [ (.r + 1 - sin 2 7)"* rf 7 (85c) 

= 2(3 4- D^j-^:) (85d) 

where the dummy variables are related by /3 = a/2 and 7 = tt/2 — a/2 
and where K and £ are the complete elliptic integrals of the first and 
second kind respectively. The argument of the elliptic integrals is the 
parameter m of Abramowitz and Stegun. 14 The parameter m is equal 
to the parameter k 2 of Jahnke and Emde or Groebner and Hofreiter. 15 ' 10 

a. 2 Differential Equations and the Power Series Expansion of U and V 
From the definitions (84) and (85) 

^ = \V. (8G) 

The differential equations obeyed by U and V are 

[(.r 2 4-.r)£, + |]r(.r) =0 (87a) 



and 

r ,r- ri il 

V(x) = 0. (87b) 



[ 



.," + ,)£, + (2, + !)£+! 



The U differential equation is verified by substituting in the defining 
relation (84a) and then reducing the resulting equation to the indentity 

« f* d sin a , , 00 x 

= / — - — — r- ttt da (88a) 

J da [x + |(1 — cos a)]' 

f T ^ cos a + I cos a - I - f cos 2 a , , oou . 

= I -. — : — 77; ttj aa. [poo) 

Jo [x + i(l - cos a)] 5 
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The 7 equation is then easily obtained by differentiation of the U equa- 
tion. 

The roots of the indicial equations of these equations are separated 
by 1 [they are and 1 in equation (87a) and —1 and in equation 
(87b)] so that the series expansion of U is of the form 



U(x) = £ UiX* where U, = U', + U}' \ In 



and 



V(x) = E Vft where 7, - 7J + V'/ \ In 



16 

x 

16 

x 



(89a, b) 



(90a, b) 



The form of the logarithmic terms has been chosen with some foresight. 
Substitution of the expansions into the differential equations and 
comparing coefficients gives the recursion relations 



jjn _ _il lL Tjr, 



» 2: 1, 



jii = J - 



i0" + i) L° 2) ' iO' + D 'J' 



J 2 1. 



»;-. -f?44)Vr 



i + i 



i ^ 0, 



(91a) 
(91b) 
(92a) 



K 1 + 1 /"» 



i + i 



0" + i) 2 _ 



' ^ 2/ ' 2 ; + 1 



7J' 



i ^ 0. (92b) 



The starting values of 7< , 7J' , ?7$ and E/J' are determined directly 
from the integral definitions of the functions (84) and (85) and the 
differential equation (86) relating U and 7. This is quite straightfor- 
ward except for 7 which must be expanded 



F « - 2 /.' f+W + 2 /.' 



and evaluated as a limit 



lim V(x) = In 



jgL + Otf,^ 



1() 



(93a) 



(93b) 



The limiting value of 7 may also be obtained quite easily from equation 
(85d) and the tabulated properties of the complete elliptic integral of 
the first kind. 15 The expansions of U and 7 are thus 

tf(*) = (2 + !* + !*'- (I** + 2fe *'+-) 
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x 
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(94) 



and 



21 



185 



r(. c ) = ^o + -,---,- + 768 



In 



ill 

x 



(95) 



A.3 Expansion of U and V in Terms of Inverse Powers 

Making a Taylor series expansion of equations (84b) and (85b) 
respectively and integrating yields the expansions of U and V in terms 
of the inverse powers of the argument 

U(x) - 2x i f' (1 + x -1 sin 2 /J)» dp 



(2j)! - (i - t) 



: J§2;- 1(-4)W 

v - 1 rc'/)! 



sin 2 ' |8 dp 



.-(»-*) 



and 



ft(2j- D(-ifi)' L (JO". • 
l'(.r) = 2x-> f (1 + i _1 sin 2 0)'* d/3 

-0 



(96a) 
(96b) 



Sfi! 

LO'O* J 



.-o + i> 



""^(-16)' 



= 7r.r 



(97a) 
(97b) 



A.4 Definition of the L„ Functions 

The L„ functions are defined in expression (34) by 



f T ( 1 — cos na) da 

lM = J U + Ki -^os"^' 



.r ^ 0, w ^ 



(98a) 
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or with the change of variable /3 = a/2 

^^fV+Sfo* ->0,«>0. (98b) 

It can be seen directly from equation (98b) that for a fixed value of n 

^^ < 0, L„(oo) = 0, and L n (x) = 0. (99a, b, c) 
ax 

From definitions (84b), (85b), and (98b) 

L,(x) = 2[U(x) - xV(x)]. (100) 

The higher L functions are determined by means of a recursion relation. 

a. 5 The L n Recursion Relation 

The L n recursion relation is 

L n+l (x) = ^r\ [4n(2.T + l)L n (x) - (2n - l)L n _,(a-) - 8nxV(x)] t 

n fc 1. (101) 

The recursion relation is verified by substituting in the definitions of 
L„ and V, equations (98a) and (85a), and reducing the resulting equation 
to the identity 

= f T 4- {smna[x + h(l - cos a)]*} da (102a) 

Jo aa 



f 



n cos na[x + h(l — cos a)] + f sin na sm a . , imu 

P — ; — rr- rTi aa. (lUJb) 

[x + |(1 - COS a)J J 

The initial functions L (a;) and L t (x) are given by equations (99c) and 
(100). Note that for large values of x the recursion relation is unstable 
for increasing n. 

A.6 Power Series Expansion of the L n Function 

The function L (x) is identically zero, equation (99c). The series 
expansion for L t (x), obtained from equations (100), (94), and (95), is 

r / \ (a jl 13 2 , 9 3 5255 , 

L,(.r) = (4 + .r - - x + ^ x - j^ * + - 



(f| _ 2x . + |. r2 _l|. c 3 + 17|,, + ...)l ln 



Hi 



(103) 
From the recursion relation (101) and the initial functions (99c) and 
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(100) the general form of L n (x) is 

L n (x) = u n (x)U(x) + v H (x)V(x) (104) 

where u n (x) and v„(x) are polynominals of order n or less in x. Because 

of the form of U(x) and V(x), expressions (89) and (90), an expansion 
of the form 

L n (x) = EL n „,r' (105a) 



where 



L Hti = U.i + L'^ \ In 



16 
x 



(105b) 



clearly exists. 

Expressions for either the coefficients in the polynominals u„(x) and 
v„(x) or the L„ ti may be determined in closed form by similar methods. 
It has, however, proven more useful to use the recursion relation di- 
rectly when the complete expression of the form of expression (104) 
is desired and the expansion (105) when a power series is desired. 

To obtain the L n ,, the expansion (105) is substituted in the recursion 
relation (101) and coefficients of x are compared to obtain a hierarchy, 
in j, of recursion relations, each member of the hierarchy being factor- 
able and depending only on the preceding member. These recursion 
relations are then factored and successively summed. 

The coefficient of x' is 

L - +1 -" = 2^T1 {4nLni ~ (2?? ~ 1)L "- 1 -' + 8n[z "-'-' " 7 '- ,]} ' 

n ^ (106a) 

where 

/,„,_, = and !'_, = 0. (106b, c) 

With the definition 

Q nti s (2n - 1)[L„., - L B _ 1>f ] (107) 

the second order recursion relation (106) factors into two first order 
recursion relations 

Q n+I ., - Q-.i + S»(L n .,_, - 7,-0, n ^ (lOSa) 
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and 

L B+l>J = L.., + 2^Vl 0n+in W " l - (1 ° 8b) 

The recursion relations for j = and j = 1 will now be summed. 
From expressions (99c) and (105a) 

Lo„ - (109) 

so that using expression (107) 

Gi.|-Iri.|. ( 11Q ) 

For j = using expressions (106b) and (106c), the recursion relation 
(108b) becomes simply 

e. + i. a = q-.o. (in) 

By inspection of expression (103) the initial value of Q„. Q is 

Qi.. = Li.. - 4 (112) 

so that from expression (111) 

Q n , = 4, n £ 1. (113) 

The recursion relation (108b) thus becomes 

L„, = L„_ li0 + gj^-j, n ^ (114) 

which with the initial value of expression (109) may be summed to yield 

°' n = °' (115a) 

4^^ n>0. (H5b) 

From the form of the expansion (105) it can be seen that 

L„(0) = L„. , (116) 

so that with (99a, b) 

-4 £ ^4-r ^ L n (.r) - L„(0) ^ 0, « ^ 1. (117) 

For evaluating the Q n>1 and L„,i sums, two relations are needed: 
2-j 2s or. _ i ~ Z-( 2^ 97- _ i 

V n — k — - 
"ft? 2fc - 1 
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= J(2« + l)(E^)-| 



and similarly 



tit^T-i^+D'S^ri 



— -n . 
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(118) 

(119) 



Evaluation of the sums of expressions (118) and (119) is analogous to 
integrating x" log x where n = and 1. With sufficient patience the sums 
can clearly be carried out for any finite n. 

For j = 1, expression (108a) becomes [using expressions (95), (109), 
and (115)] 



Q.+ ..1 = Q..i +8n 



- 4 S^ 



- In 



16 



n > 1 



(120) 



where the initial function is [using expressions (103), (105), and (110)] 



Q U1 = L 1A = 1 - In 
Summing [using expression (119)] yields 
Q. + i.i = 1 - In 



Hi 



(121) 



^ + E 8* 

# i • - 1 



- -(2n+l) 2 ln 



16 

x 



.4- 1 i I 1G II 

+ 4(^+1)- Eg^ -(4^-1), 



The L„.i recursion relation is then 

Lw... - />,.., + (2n+l)[- In' 1 " 



n £ 1. (122) 
n £ 1 (123) 



which may be summed using the initial value of expression (121) and 
the sums (118) and (119): 



&»+,., = 1 - In 



- I [l + E <fr + 1)1 






(124) 
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&... = 



0, n = 


0. 


-2«°(i in 


X 



)_ 2 ^ + (4n »_ 1) |: 2T i- r , 



(125a) 



n ^ 1 . (125b) 

It is clear that the procedure leading to expressions (115) and (125) 
may be carried onward to lead to an expansion of the form 



oo r 



16 



L n (x) - L„(0) = + £ £ L (1) (;, k) § hi ±2 + L «>(j t k) 



+ 



L-<f,« t^H 



(126) 



where the L <p> 0", k) are functions of ; and A; only. It is clear that in the 
limit x — > 0, n — » co an expansion in terms of 



7T 

:m 2 



(127) 



may be made where \/h is introduced as the finite expansion parameter. 
Replacing the sum in expression (125) by its approximating integral 
yields 



where 



lim [L n (x) - L n (0)] = -2tt 2 ~ In 

._« A 



n = irhx *X~\ 



4eX 

7r/l 



+ 4 (^) (128a) 



(128b) 



A.7 Expansions of L m Terms of Inverse Powers 

The expansion of L in terms of inverse powers of the argument is 
obtained by Taylor expansion of expression (98b) and integrating. This 
yields 

L ^ )=4X i (l+^sin^^ 



= 4Z 



(2j)! 



jV.M\2 ^ 



ft (-4 )-iy!i 
where for example 






t/2 



sin 2 n/3 sin 2 ' /3 e//3, 



£.(*) - 



f, (2j + 1) 1 r(2j)!l 2 T -" + »> 
* &! 0+D (-16)'L0!) 2 J 2 



n ^ 1 (129a) 
(129b) 
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By considering the Fourier decomposition of sin' /3 it can be seen that 

/»t/2 , nr/2 

I sin 2 n$ sin 2 ' d& = r / sin 2 ' d/3, n > ; 

Jo 4 

(independent of n) so that 

£o(-c) - 0, 

*.<»> = » £ pfe? Ca&G"*" < '** > + 0[x ~ <n+i}]) ' n > [ (131a) 

or identifying with expression (97) 

L n (x) = V(x) + 0[*" <m+|, l, » > I- (131b) 

Evaluating expression (129a) directly in those cases in which expressions 
(129b) or (131a) cannot be used yields 

L Q (x) = 0, (132a) 

L.Cr) = rx- u * - ^ x~ m + ^ x- B/2 + 0(.r- 7/2 ) (132b) 

L 2 (.r) - «-" - \ x-*' a + H aT" + OCT""), (132c) 

L n (.r) = «-"■ - | af 8/a + || x" 5/2 + 0(af 7/2 ). (132d) 

a.8 The Gaussian Transformation 

In the neighborhood of £ = 1, the convergence of the power series 
in a: or aT 1 is rather slow. Either the gaussian or Landen transformations 
may be used to transform the U, V, or L x functions into a region of rapid 
convergence. 10 In the present case, the gaussian transformation is pre- 
ferred since it does not introduce incomplete elliptic integrals as does 
the Landen transformation. 

The result of the gaussian transformation is 

x, = 4x»(l + *)*[(1 + x? + x*] 2 , (133a) 

or inversely 

x\ 
X = 4(1 + *,)*[(! + *) f + I] 2 (133b) 
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for the argument and 

V(x) = 2TV(x 1 ), 

TU(x) = U(x t ) - | 7(.r,), 

TL t (x) = InW + j-2^5 V(x 1 ), 

for the functions where 

T = (1 + xrf = (1 + a;)* + xK 



(134) 
(135) 

(136) 
(137) 



APPENDIX B 

The Force and Stability Functions 

This appendix is a compilation of expressions for the force F and 
stability S n functions. Each of the functions is written in terms of the 
U and V or L n functions of Appendix A. Expressions in terms of the 
complete elliptic integrals of the first and second kind (denoted by K 
and E respectively) permit the use of tables 1 ' 1 or numerical computation 
using the Landen transformation or the gaussian transformation. 16 
The gaussian transformation is used in Section A.8. The power series 
expansions provided are necessary in obtaining numerical values of the 
functions for either very large or very small values of the argument and 
also provide the asymptotic forms of the functions. The argument of 
the functions is the domain diameter to thickness ratio, d/h = 2r /h 

b.i The Force Function 

The force function is written in terms of U by comparing the form 
of F, expression (40), and the form of U, expression (84b), 



F Xh) r \h 



U 



:- _ 9 



(138a) 



This expression is written in terms of the complete elliptic integral of 
the second kind using expression (84d) 



F (i) = 2 (i 

W t \h 



i + j£Yb[(i + h 2 /d 2 y l ] 



■]■ 



(138b) 



It is expanded about h/d = using expression (94) 

665 



it) - 1 



2 + 32 \d) 64 \dl + 24576 \dl "* 
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In Uf 

h 



(d\ m d _ 2 Id 

r \h h r \h 



(138c) 

Additional terms may be generated using expressions (89) and (91). 
It is expanded about d/h = using expression (96b) 

Additional terms may be generated using expression (96a). 

B.2 The Radial Stability Function 

The radial stability function is written in terms of U and V using the 
definition of the radial stability function [equation (66a)], the expression 
for F [equation (138a)], and comparing the derivative of F [equation 
(41)] with the form of V [equation (85b)], 



<f)--l(i 
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or in terms of L x using the expression for L x [equation (100)] and the 
expressions for L x (0) [equations (115) and (116)] 



S °{h) ~ 2tt U 



Mf J - U0) 



(139b) 



Expression (139a) is written in terms of the complete elliptic integrals 
of the first and second kind using expressions (84d) and (85d), 
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The expansion about h/d = is obtained using expressions (139b) and 
(103). 
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Additional terms may be generated using expressions (139a) and (89) 
through (92). The expansion about d/h = is obtained using expressions 
(139b), (132b), and (103) to obtain L, (0) 
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(139e) 



Additional terms may be generated using expression (129b). 

b.3 The Elliptical Stability Function 

From the general definition of the S n of expression (66b) the elliptical 
stability function is 



«© - "I 



*©+£(ffM*i -'■<■»> 



(140a) 



Using expression (139b) for £ > the L n recursion relation (101) to reduce 
L 2 to Lj , and V and (103) to obtain Z»i(0), *S 2 is written in terms of I/! 
and 7: 
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The function L x is then eliminated using expression (100) to obtain the 
expression in terms of U and V: 

(140c) 

which then is written in terms of the complete elliptic integrals of the 
first and second kind using expression (84d) and (85d) : 
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The expansion about h/d = is obtained using expressions (140b), (95), 
and (103): 
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Additional terms may be generated using expressions (140c) and (89) 
through (92). The expansion about d/h = is obtained using expressions 
(140b), (97b), and (132b): 

*«)-*«)■-£«)■+£«)'+■•■■ ^ 

Additional terms may be generated using expressions (140b), (97a), and 
(129b). 

B.4 The General Stability Functions 

Using the definition (66b) and the expression for S [equation (139b)], 
the S n are written in terms of the L n as 



a (i\ - i jl ($? 

"W ' n 2 - 1 2tt \h 
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n ^ 2. (141a) 



The leading term of the expansion about h/d = is obtained using ex- 
pressions (103), (105), (115), (116), (125), and (126): 

(d\ 1 f. 4d An 2 - 1 ^ i_ 2 n 2 + 1 1 /w/A 
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The expansion about d/h = is obtained using expressions (115), 
(116), (129b), and (131a): 
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APPENDIX C 

Symbol List 

Numbers in parentheses are defining equations or figures. 

A exchange constant (71) 

a area 

d mean domain diameter, 2r (42) 

E(x) complete elliptic integral of the second kind (x = m = k 2 ) 

E H energy due to applied field (9) 

E M internal magnetostatic energy (10) 
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E T total energy (7) 

E w total wall energy (8) 

F(x) generalized radial force function (33, 138, Fig. 3) 

H magnetic field vector 

H uniform applied field 

H a anisotropy field 

H N nucleation field (74) 

(H,) RV z-averaged z-component of magnetic field (45) 

h plate thickness (Figs. 1, 4) 

K(x) complete elliptic integral of the first kind (x = m = k 2 ) 

K u uniaxial anisotropy constant (71) 

L n (x) integrally defined function (34, 98) 

I characteristic length, ov/4-n-M* (67) 

l„ wall width (81) 

M magnetization vector 

M, saturation magnetization 

n rotational periodicity (1) 

k terms of order k 

q quality factor, K u /2tM 2 , = HJ4wM. (72) 

r cylindrical coordinate (Figs. 1, 2) 

r, plate radius (Fig. 3) 

r n nth radial Fourier amplitude (1) 

r„ mean domain radius (1) 

S n (x) nth infinite plate stability function (66, 139, 140, 141) 

s distance between interacting magnetic charges (10b, 24, 28) 

U(x) integrally defined function (84) 

u(x) unit step function (5) 

V volume 

V(x) integrally defined function (85) 

z cylindrical coordinate (Figs. 1, 4, 5) 

Z operator (23) 

z z-z' (21) 

A-E variation in energy (11, 43) 

Ar„ variation in r„ (3) 

A0 n variation in 0„ (3) 

8(x) dirac delta function 

f d' - 6 (27, Fig. 2) 

7] polar azimuthal angle (Fig. 5) 

■q w polar angle of wall normal (Fig. 5) 

6 cylindrical coordinate (Fig. 1) 

d n nth Fourier phase angle (1) 
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X n wavelength of nth variation (54) 

v azimuthal angle (Fig. 5) 

i> u azimuthal angle of wall normal (Fig. 5) 

£o wall displacement vector (Fig. 5) 

p coordinate in displaced cylindrical coordinate system (35, Fig. 2) 

<r„ wall energy density 

cp coordinate in displaced cylindrical coordinate system (35, Fig. 2) 

Xt transverse susceptibility 

£2 magnetostatic potential 
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